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Abstract. - We study in detail by means of ab-initio pseudopotential calculations the elec- 
tronic structure of five- fold coordinated (T5) defects in a-Si and a-Si:H, also during their for- 
mation and their evolution upon hydrogenation. The atom-projected densities of states (DOS) 
and an accurate analysis of the valence charge distribution clearly indicate the fundamental 
contribution of Tg defects in originating gap states through their nearest neighbors. The 
interaction with hydrogen can reduce the DOS in the gap annihilating T5 defects. 



The atomic origin of the midgap energy levels in a-Si has been commonly ascribed to 
three- fold (T 3 ) coordinated atoms, and the observed reduction upon hydrogenation has been 
Q , explained with the passivation of the "dangling bonds" by H jjj . This picture has been pro- 

posed in analogy with the well known mechanism occurring in undercoordinated configurations 
such as at surfaces and vacancies in crystalline silicon (c-Si) g. 

The dangling bond picture for a-Si has been widely supported by many electronic structure 
calculations || [|, g, |], a- although it has been recognized that gap states can be induced 
also by other coordination defects (j, f?j, e.g. five- fold (T 5 ) coordinated "floating bonds" and 
anomalous four- fold (T4) coordinated atoms. 

We focus our attention here on the role of T$ defects. Their importance in a-Si has been 
clearly stated by Pantelides |g, || and Kelires and Tersoff JH]] a douzen of years ago. Pantelides 
[|| || suggested that the presence and the role of T 5 defects must be seriously reconsidered in 
order to explain some theoretical and experimental data (effective electron correlation [[ll|, 
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hyperfine structure from electron paramagnetic resonance data [{l2[, relationship between 



intensity of paramagnetic signal and density of Si-H bonds 13 ) that otherwise would remain 
unexplained in the common picture involving only dangling bonds, and he gave some arguments 
suggesting that T 5 defects could be predominant. He argued that T 3 and T 5 are conjugated 
defects, since a bond elongation can transform a T 4 + T 5 structure into a T 3 + T4 one B; 
furthermore, he proposed a mechanism for H diffusion based on floating-bond switching and 
annihilation/formation of T 5 's through interaction with H M. The empirical simulation by 
Kelires and Tersoff [[10| has shown that T5 atoms have lower energy than T3 atoms, and 
therefore should be favoured in general. Also some ab-initio molecular dynamics simulations 
of a-Si structures indicate a predominance of T5 defects with respect to T3 H, |15[ . 

The ideas of Pantelides have then been applied mainly in discussing the geometrical char- 
acterization of defects || 0, [Uj [IT], [y|. Our aim is to discuss their soundness in terms of 
electronic properties, by analyzing charge distributions and atom-projected DOS obtained 
from accurate ab initio calculations, also following some possible process of formation of T5 
defects and of evolution upon hydrogenation. 

To this purpose, we start from some selected samples generated by other authors []14[ [l5, 19 
using Car-Parrinello molecular dynamics (CPMD). These structures are a good starting point 
for this study, since they contain floating bonds; furthermore, they reproduce quite well the 
experimental pair correlation function and bond angle distribution function using a reasonable 
number of atoms and hence they are suitable for accurate ab-initio studies. The small number 
of defects, which is however larger than in experiment, allows to easily single out effects 
associated to local features. The configurations studied are cubic supercells of side a = 2 a , 
where ao is the equilibrium lattice parameter of c-Si. With respect to the original configuration, 
where ao was fixed to the experimental value, we use the theoretical equilibrium lattice 
parameter ao — 10. 17a. u., which also corresponds — in our calculations — to the optimized 
density of a-Si and a-Si:H. The starting configurations contain respectively 64 Si atoms to 
describe a-Si @, [l5| and 64 Si atoms plus 8 H atoms for a-Si:H @ |§. 

We have studied both the mean configuration at room temperature and a snapshot of 
the CPMD run, in order to check for possible anomalies due to statistical average, and also 
to single out the effect on the electronic features of tiny structural variations. The CPMD 
configurations, aiming mainly at reproducing the structural properties, have been obtained 
using a kinetic energy cutoff E cut —12 Ry and the T point only for Brillouin Zone (BZ) sampling 
[fbjl [19, H3]. Since electronic structure studies require better accuracy, we improve in our 
calculations the BZ sampling using 4 inequivalent special k points for self-consistency and 75 
k points for DOS. These parameters have been chosen as a reasonable compromise between 
accuracy and computational cost, after tests with E cut =16 Ry and with the T point or 32 
special k points for self-consistency. We have used for Si the pseudopotential by Gonze et al. 
[ p0| and for H a smoothed Coulomb potential. 

The optimization of the a-Si and a-Si:H structures with the new computational parameters 
is accompanied only by small structural rearrangements. The results for the structural and 
electronic properties of the mean relaxed configurations that we present here are essentially 
valid also for the others (snapshot, unrelaxed). 

The mean structural properties of these configurations are similar to those discussed in refs. 
[[ill MS, pj| . We only report here that in a-Si the mean bond length is d ~ 4.47 a.u., quite 
similar to the crystalline one which is 4.40 a.u.. The location of the first minimum of the 
radial distribution function defines geometrically the cutoff distance for the nearest neighbors 
(NN), which for Si-Si turns out to be Rnn = 5.08 a.u. and Rnn = 5.49 a.u. in a-Si and 
a-Si:H respectively. In a-Si:H each H is bound to one Si atom with an average distance du = 
2.95 a.u.. With those values of Rnn, the resulting average coordination number for Si in a-Si 
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Fig. 1. - Total and projected DOS of c-Si, a-Si and a-Si:H. DOS is normalized to unity and 
projected DOS to the number of valence electrons for atomic type. The energy scale is referred 
to the Fermi level. Upper panel: total DOS for a-Si (solid gray line), a-Si:H (solid black line), 
and c-Si (dashed line). Middle panel: atom- projected DOS for a-Si on T5 atoms (dashed line), 
T4 atoms which are NNs of T5 (solid black line) , and the remaining T4 atoms (solid gray line) . 
Lower panel: atom-projected partial DOS for a-Si:H on T4 atoms which are NNs of T5 (solid 
black line), T5 atoms which are NNs of T5 (dashed line), the remaining T4 and T5 atoms (solid 
gray line), hydrogen atoms (dotted-dashed line). 

Fig. 2. - Partial atomic charge corresponding to the occupied states in the gap of a-Si. Almost 
all the relevant contributions are associated to some NNs of the T5 atoms, i.e. to atoms labeled 
4, 11, 12, 35, 40, 49 and 59, which are NNs of atom 13 or 53, or of both, as it is the case of 
atom / = 40. The peculiar local structure is sketched. "Long" bonds (see text) are in black. 



and a-Si:H is slightly larger than 4. For the purpose of the present work a classification of 
coordination defects based on a pure geometrical analysis is sufficient, but we point out that 
possible ambiguities in the bonding configuration can be resolved only by an accurate analysis 
of the valence charge distribution and of the "electron localization function" [Ell E2| . 

Both our a-Si and a-Si:H configurations have a metallic character, with coordination-defect 
induced gap states (the calculated total DOS is reported in the upper panel of fig. 1). The 
precise definition of the energy gap is one of the unresolved problem for an amorphous system. 
We have chosen to identify bottom (Eq*) and top (E^ p ) edges of the gap with the absolute 
minima of the DOS below and above the Fermi level Ep. The occupied states in the gap are 
the most relevant for our discussion, even if their contribution to the charge density is small 
(less than 1 % for a-Si and 0.5% for a-Si:H). Their atomic origin can be identified through 
the projection on LCAO pseudo-wavefunctions. The accuracy of the projection, evaluated by 
means of the spilling parameter S |2j| measuring the incompleteness of the LCAO basis set, 
is satisfactory (S'=0.03). 

We focus now our attention on a-Si. The configuration considered here contains two T§ 
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sites close one to each other, with an "interstitial" (J) atom connecting them (see fig. 2), and 
no "true" T3 defects. A charge density analysis helps in characterizing the different types of 
bonds [[15], gl[ . We observe that T$ sites are accompanied by a valence charge density depletion. 
In particular, those T5-T4 bonds which are about 10% longer than the average value ("long" 
bonds) and the bonds T 5 -I are characterized by a charge density distribution remarkably 
smaller (about 1/3) than the corresponding one in the crystalline phase, or than the average 
in the amorphous phase. These bonds are therefore weak and are the best candidates to break 
under a network distortion, giving raise to a T4 + T§ — > T3 + T4 transformation. 

The middle panel of fig. 1 shows the DOS projected on different classes of atoms: T5 atoms, 
atoms which are NN of T5, and the remaining T4 atoms. The peak of the DOS in the gap is 
unambiguously associated to the NNs of T5 sites. In fig. 2 we show the atomic projection of the 
charge density of the gap states, obtained again from the LCAO coefficients. It is even more 
clear that the gap states are not located on T^s themselves but are delocalized on their NNs, 
and in particular on those connected by longer and weaker bonds. Evidence of a midgap state 
with these characteristics is also reported in ref. |l8| . We argue that previous works ascribed 
the origin of gap states uniquely |3| or mainly H^f^, |23|, |24j to T 3 , simply because they 
disregarded an accurate analysis of NN environments. 

From inspection of fig. 2 we observe the large contribution to the DOS in the gap coming 
from the atom / which is NN of two T$'s. We also observe that also some normally coordinated 
sites (labelled 18 and 44), characterized by having "long" bonds, partially contribute to gap 
states, in agreement with similar findings reported in refs. [0, p], B3 |. 

The a-Si:H |19|] structure studied here is more complex. Well defined T5 sites are present 
and are the predominant kind of defect (we do not address in this paper the discussion of 
other more complicated defects which could require a better characterization beyond the pure 



geometrical analysis |2lj, 22 ). Overcoordination arises from five Si atoms, or from four Si 
atoms and one hydrogen as NN. The main features of the electronic structure are not very 
different from those of a-Si. As we found in a-Si, those atoms which are NNs of some T5 sites 
give a large contribution to the DOS in the gap (fig. 1, lower panel). 

In order to get further insight onto the characterization of the electronic properties of T5 
defects and to recover the ideas of Pantelides, we want to follow in terms of electronic states 
some possible processes of (a) T3 — > T5 transformation and (b) annihilation of T5 by H. To this 
purpose, we follow two "computational experiments" starting from the CPMD configurations 
described before. 

In the first one, we start from a-Si:H sample. We focus our attention on a normally 
coordinated cluster of atoms (fig. 3, upper panel) containing two T4 Si (labelled 1 and 4) 
close one to each other — but not directly connected — and their NNs which include two H 
atoms. This structure T4 + T4 has an almost negligible DOS in the gap. We first create 
artificially a T3 site by removing one H atom and leaving all the other atoms fixed, so that 
the structure is now T4 + T3 (fig. 3, central panel). This configuration corresponds to a huge 
DOS in the gap, almost entirely originating from the dangling bond. After relaxation of the 
atomic positions the distance between atoms 1 and 4 is reduced from 5.63 a.u. to 4.48 a.u. 
and a new bond is formed, so that the structure is now T5 + T4 (fig. 3, lower panel). This 
process (removal of one H followed by bond reconstruction) is similar to that studied in ref. 
[181, with the difference that the other H atom in our configuration is replaced there by one 
Si atom. We observe that in the final configuration (T5 + T4) there are still midgap states, 
but now they are delocalized on the NNs of the atom 4 and pushed towards lower energies. 
Our findings support the picture of conjugated T3 and T5 defects. Moreover, the spontaneous 
evolution from T4 + T3 to T5 + T4, which is in agreement with the findings of ref. |lq ], confirms 
-although in a particular case — that the T$'s can be energetically favoured with respect to 
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Fig. 3. - Projected DOS of a cluster of Si (gray) and H (white) atoms selected from a-Si:H, 
evolving from a T4 + T4 (upper panel) to a T4 + T3 configuration (central panel) after removal 
of one H atom, and finally to a T5 + T4 configuration (lower panel) after optimization of atomic 
positions. Solid black line is for atom 1, dashed line for atoms 4, 13, 16, 58 (whose contribution 
to the DOS remains essentially unchanged), grey line for atoms 59, 10, 50, 65, which are NNs 
of atom 4. 

Fig. 4. - Projected DOS of the cluster of atoms selected from a-Si initially containing the two 
T5 defects (upper panel) , and after the insertion of a first H nearby the T5 labelled 53 (central 
panel) and a second one nearby the T5 labelled 13 (lower panel), after optimization of atomic 
positions. Dotted line is for atom 40, solid black line is for the remaining NNs of atom 13, 
dashed line for NNs of atom 53, thin gray line for atoms 13 and 53, whose contribution to 
the DOS in the gap is negligible and remains essentially unchanged in the evolution of the 
structure. 



the T 3 's, as indicated in ref. [flCfl - 

The inverse process, that may occur in real systems as proposed in ref. M, indicates a 
possible mechanism of annihilation of T5 defects by H. But this can be better studied following 
in terms of electronic states our second computational experiment. We start from the a-Si 
configuration and we focus our attention on the cluster with the two T5 sites (fig. 4, upper 
panel). We add one H atom nearby a T5. Relaxing the whole structure, a T4 — T5 "long" bond 
breaks, so that a new Si-H bond is formed and the T5 becomes a normal T4 site (fig. 4, middle 
panel). This is accompanied by a lowering of the DOS in the gap. Adding another H atom, 
also the second T5 is annihilated and the final structure has a vanishing DOS in the gap (fig. 4, 
lower panel). 

A variety of other structural models have been proposed in the literature for a-Si and 
a-Si:H, describing even better some global properties, but we are confident that our findings 
concerning local environments are not much affected by the particular choice of the structural 
model. Infact preliminar results of another work in progress p2| using a larger sample of a-Si:H 
generated by L. Colombo using tight-binding molecular dynamics pq confirm the validity of 
the present findings. 
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In conclusion, we have presented the results of accurate ab initio self-consistent pseudopo- 
tential calculations of the electronic properties of T5 defects in a-Si and a-Si:H starting from 
some configurations generated via ab initio molecular dynamics, and we have followed in terms 
of electronic density of states some possible processes of their formation and annihilation by H. 
Our calculations clearly indicate that i) T5 sites give a fundamental contribution to gap states 
through their NNs (fig. 1,2); ii) in some circumstances the evolution from T3 to T5 is favoured, 
and this is accompanied by a derealization of the corresponding electronic states with a shift 
towards lower energies (fig. 3); Hi) the DOS in the gap can be reduced by the interaction of H 
with the T 5 's (fig. 4). These results support Pantclidcs' proposals, and suggest that T 3 and T 5 
in a-Si not only are "structurally" conjugated, but they also play a similar role in originating 
gap states, which can be in both cases passivated by hydrogen. This pictures reconcile those 
experimental findings that a description uniquely emphasizing the dangling bonds would not 
explain. 
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(SdG) acknowledges support from the MURST within the initiative "Progetti di ricerca di 
rilevante interesse nazionale". We would like to thanks L. Colombo and G. L. Chiarotti 
for useful discussions and suggestions. 

REFERENCES 

[1] Davis E.A., Journal of Non-Cryst. Sol., 198-200 (1996) 1; Ley L., in The Physics of Hydro- 
genated Amorphous Silicon II, edited by Joannopoulos and Lucovsky, (Springer) 1984 p. 61. 
[2] Pickett W.E., Phys. Rev. B, 23 (1981) 6603. 
[3] Fedders P. A. and Carlsson A.E., Phys. Rev. B, 37 (1988) 8506; Phys. Rev. B, 39 (1989) 

1134. 
[4] Biswas R., Wang C, Chan C, Ho K. and Soukoulis C, Phys. Rev. Lett, 63 (1989) 1491. 
[5] Holender J., Morgan G. and Jones R., Phys. Rev. B, 47 (1993) 3991. 
[6] Lee In-Ho and Chang K.J., Phys. Rev. B, 50 (1994) 18083. 
[7] Tuttle B. and Adams J.B., Phys. Rev. B, 53 (1996) 16265. 
[8] Pantelides S.T., Phys. Rev. Lett, 57 (1986) 2979. 
[9] Pantelides S.T., Phys. Rev. Lett, 58 (1987) 1344. 
[10] Kelires P.C. and Tersoff J., Phys. Rev. Lett, 61 (1988) 562. 
[11] Bar- Yam Y. and Joannopoulos J.D., J. Electron. Mater., 14 (1985) 261. 
[12] Biegelsen D.K. and Stutzmann M., Phys. Rev. B, 33 (1986) 3006; Stathis J.H., Phys. Rev. 

B, 40 (1989) 1232. 
[13] Biegelsen D.K., Street R.A., Tsai C.C. and Knights J.C., Phys. Rev. B, 20 (1979) 4839; 
Kumeda M. and Shimizu T., Jpn. J. Appl. Phys., 19 (1980) L197. 

[14] Buda F., Chiarotti G.L., Stich L, Car R. and Parrinello M., J. Non-Cryst. Sol, 114 
(1989) 7. 

[15] Stich L, Car R. and Parrinello M., Phys. Rev. B, 44 (1991) 11092. 

[16] Fedders P.A. and Carlsson A.E., Phys. Rev. Lett., 58 (1987) 1156. 

[17] Fedders P.A, Drabold D.A. and Klemm S., Phys. Rev. B, 45 (1992) 4048. 

[18] Tuttle B. and Adams J.B., Phys. Rev. B, 57 (1998) 12859. 

[19] Buda F., Chiarotti G.L., Car R. and Parrinello M., Phys. Rev. B, 44 (1991) 5908. 

[20] Gonze X., Stumpf R. and Scheffler M., Phys. Rev. B, 44 (1991) 8503. 

[21] Peressi M., Fornari M., de Gironcoli S., De Santis L. and Baldereschi A., Phil. Mag. 

B, in press. 

[22] Fornari M., Peressi M., de Gironcoli S. and Baldereschi A., unpublished. 



M. FORNARI et al.\ FLOATING BONDS AND GAP STATES IN a-SI ETC. 7 

[23] Drabold D., Fedders P., Sankey O. and Dow J., Phys. Rev. B, 42 (1990) 5135. 
[24] Holender J. and MORGAN G., J. Phys.: Condensed Matter, 4 (1992) 4473. 
[25] Sanchez-Portal D., Artacho E. and Soler M., J. Phys.: Condensed Matter, 8 (1996) 3859. 
[26] Colombo L., in Ann. Rev. Comp. Phys., IV, edited by Stauffer D., (World Scientific, Singa- 
pore) 1996 p. 147. 



